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Abstract 

A constitutive model for the stress-strain relationship of single forest root systenn was developed in order to provide 
theoretical foundations for the mechanisnns of soil-reinforcement by root system and offer a reliable basis for the analysis of 
root tensile strength character. This study started a general form of linear and non-linear stress-strain relation that was 
mathematically defined by four boundary conditions observed in typical tensile tests of single roots. The parameters of the 
model were determined by experiment data and had definite physical meaning. The model was verified by experiment 
data, which showed that the calculated values were in good agreement with the experimental single root tensile test 
results. The constitutive model was validated and found to be feasible for modeling single root tensile stress. 
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Introduction 

Roots are liey plant organs tliat perform vital functions such as 
absorbing and storing soil moisture, nutrients, and energy, as well 
as anchoring the plant to the soil. One of the main ecological 
functions of forest roots is soil stabilization, which can enhance soil 
conservation, stabilize slopes, prevent gravity erosion, and prevent 
river and reservoir bank scouring [1-4]. Since 1990s, the 
mechanism of soil-reinforcement by root systems has received 
widespread attention from research community and led to many 
interdisciplinary studies [5-7] . Since single root of plants is the 
basic unit of soil-reinforcement, a series of studies has been carried 
out to study the mechanism of soil-reinforcement by single root 
system [8]. 

Currently indoor root tensile tests and in situ root-soil composite 
horizontal and the vertical puUout tests are used to determine the 
tensile strength, shear strength, and factors contributing to the root 
system's resistance to shear of different species [9-10]. Some 
models have been developed to describe the constitutive relation- 
ship(the stress-strain relationship) of single root system. Most 
models use different mathematical functions, such as linear, 
hyperbohc, polynomial, or power function [11-13] to describe the 
relationship between stress and strain. Those functional constitu- 
tive relationships are empirical in nature, since mathematical 
functions are selected based on die best curve fitting to 
experimental data. The experiment data determine the constants 
(or parameters) of the mathematical functions, which do not have 



clear physical meaning. The applicability and accuracy of those 
models are gready limited [13-14]. 

In this study, boundary conditions that a constitutive model for 
tensile stress of a single root system should conform to are 
investigated by analyzing typical stress-strain curve of single root 
systems. A new constitutive model for the tensile stress resistance of 
single root systems is proposed and mathematically solved using 
these boundary conditions. Furthermore, the term of critical point 
(the definition will be presented in Section 2.3) is introduced and 
derived from the model. Finally, the model is verified by 
comparing the model calculated values to experimental data, 
which show good agreement between two. This study may provide 
theoretical foundations for the mechanisms of soil-reinforcement 
by root system and give a reliable basis for the analyses of the root 
tensile strength. 

Materials and Methods 

2.1 Ethics Statement 

The sampling site is managed by Mulan Weichang, a state- 
owned forest Authority of He-bei province, China. Our study 
complies with the current laws of China and international rules. 
AH necessary permits have been obtained before the actual 
sampling. The field study does not involve any endangered or 
protected species. Data will be made available upon request. 
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2.2 Sample Collection and Testing 

The field samples were coUected in He-bei province, which 
surrounds the Capital City of Beijing, at 36° 01' to 42° 37' N 
latitude and 113° 31' to 119° 53' E longitude. He-bei province 
covers a total area of 187,700 square kilometers and has a 
continental monsoon climate, with cold and dry winters and hot 
and humid summers. The annual precipitation is 400 to 800 mm, 
occurring mostly in summer. Roots of four common species of tall 
trees were chosen to test the proposed model. The trees were Pinus 
tabulaefonnis, Betula platyphjlla, Quercus mongolica and Larix gmelinii. 
Tensile tests were carried out for various root diameters and root 
lengths. To preserve the root and ensure a true reflection of the 
root mechanical properties, all root samples collected from the 
field were placed in sealed bags and stored in a refrigerator at 4°C 
before the test was taking place. All root samples were tested 
within one month and no pretreatment before testing. 22888 root 
samples in four gauges (50 mm, 100 mm, 150 mm, and 200 mm) 
were tested in total; among them, 6462 samples were in gauge 
50 mm. 

TensUe test was performed on WDW-IOOE electro-universal 
tester (Time Shijin, China; see Figure 1) equipped with 
microcomputer control. The test force range was 400-100 kN, 
and the speed range was 0.001-500 mm/min. The machine has 
fiinctions of full automatic shift and stepless speed regulation. Its 
accuracy in measuring test load and displacement is ±0.5%. 

This study uses and reports experiment data of single roots with 
different diameters and frxed gauge (50 mm) only. Root samples 
with intact bark and relatively uniformed diameter were selected. 
Each root was divided into four pieces; each piece was roughly 
60 mm in length. The diameter at three dividing points was 
measured and the average diameter was recorded as the root 
diameter. The universal tester was adjusted to a unified standard 
distance of 50 mm. Root sample was inserted (50 mm at each end 
was frxed) onto the instrument and pulled (10 mm/min) until the 
root is completely puUed apart. Necking was observed during tests. 
Necking, in engineering and materials science, is a phenomenon 
that the tensile deformation (strain) has not reached to its extreme 
when the stress reaches the maximum. The sample statistics of four 
tree species are summarized in Table 1. 

The test force (stress) and displacement (strain) data were 
coUected by measuring the force and distance required to pull 




Figure 1. Photo of the universal tester used for tensile test. 
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apart a single root. The data were then plotted as stress-strain 
curves. 

2.3 Boundary Condition 

Before the discussion of the boundary condition, the dimen- 
sionless coordinates are defmed for convenience: 

Where £ is Strain; a is Stress; e.p is Peak strain and (Jp is Peak stress. 
The following terms are also used in this study: 

E = Ea/Ep, Ep = ap/Ep, Eo=da/ds\^.^Q. 

Where Eo is Linear elastic modulus; Ep is Peak secant modulus and 
E is Dimensionless elastic modulus. Since Eg is larger than Ep, E is 
larger than 1. 

A typical experiment stress-strain curve of single root is 
presented in Figure 2. Figure 2 uses dimensionless strain and 
stress for x and y axis, respectively. It shows that a stress-strain 
curve starts from zero stress and strain. As the load (stress) 
increases, the root deformation (strain) increases linearly in the 
initial stage. As the load continuously increases, the root 
deformation remain increasing, but the slope of the curve 
decreases monotonically. The curve is concave down and 
gradually plateaus to the maximum stress at which the root 
breaks. A critical point is marked on the curve in Figure 2. Before 
the critical point, the root linear deformation is dominated and 
after the critical point, the non-linear deformation becomes 
significant. 

Above observations can be mathematically expressed as: 



P=dy/dx=0 




0.4 0.6 

Strain/peak strain 



1.0 



Figure 2. Typical stress-strain curve. (1) y=0, when x=0. The stress- 
strain curve passes the origin. (2) dyldx=-^ 1 — = -^, when x=0. The 

«£ Ep Ep 

slope of the curve at the origin is equal to the dimensionless linear 
elastic modulus of the root. (3)y= 1, whenx= I.The peak stress happens 

dy 

x= 1 . (4) — = 0, when x = 1 . The slope of the curve at the peak stress is 
ax 

zero. 
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1) j)' = 0, when x = Q. The stress-strain curve passes the origin. 

2) dy/dx = I — = -pr^, when x-Q. The slope of the curve at 

d& Zp Ep 

the origin is equal to the dimensionless linear elastic modulus of 
the root. 

3) J = 1 , when x=\. The peak stress happens x = 1 . 
dy 

4) — = 0, when x=\. The slope of the curve at the peak stress is 
zero. 

These four mathematical descriptions are the very boundary 
conditions that the proposed constitutive model should conform 
to. 



2.4 Constitutive Model 

As Figure 2 indicates, the root deformation of a single root is 
composed of both linear and non-linear deformations. Equation 
( 1 ) is a general form of linear and non-linear deformation of single 
root. 



y = ax +f.\' + (5?. 



(1) 



The four parameters of Equation (1) can be determined by 
applying the boundary conditions of (1) to (4) described in 2.3. 
Equation (1) is rewritten as: 



y = Ex-{E-\)xE=\. 



(2) 



Equation (2) is our proposed constitutive model for the stress-strain 
relationship of single root. The first term on the right hand side in 
Equation (2) describes the linear portion of root deformation, 
while the second term describes the non-linear portion of root 
deformation. 

The second derivative of Equation (2) is: 



d'-y 
dx^ 



E-\ 



2-E 
XE-l. 



(3) 



Since E is larger than 1, Equation (3) is less than zero, which 
indicates the curve is concave down, as observed in Figure 2. 

As discussed in Section 2.3, the root deformation of single root 
system is mainly linear in the early stage and non-linear in the late 
stage, the critical point is the point that separates two stages. The 
proposed constitutive model of Equation (2) represents both linear 
and non-linear deformation. Strictiy speaking, there is no absolute 
linear or non-linear deformation. In the early stage, the non-linear 
deformation, second term of the right hand side of Equation (2), is 
very small when compared to the linear deformation of the first 
term. As the curve passes the critical point, the non-linear 
deformation becomes more and more significant. 

2.5 Results and Discussion 

Data of tensile tests are used to determine physical parameters 
of single roots. Eg and Ep, used in the proposed constitutive model 
of Equation (2) for four different root species. The linear elastic 
modulus (Eg) is the ratio of the material tensile stress to the 
corresponding tensile strain in the linear range and used to 
characterize the material resistance to linear deformation capacity. 
For instance, the linear elastic modulus of Pinus tabulaeformis root 
system is 50% to 70% of the tensile strength limit [15]. It is a safe 
assumption that the root deformation is mainly linear at 40% of 
the tensile strength limit for Pinus tabulaeformis root system. 
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Therefore, the secant modulus (the ratio of the stress and strain) at 
40% of the tensile stress limit was taken as the elastic modulus 
for Pirns tabulaeformis. The same method is also 
applied to determine Eg for other three root species. The peak 
scant modulus, Ep = (jp/Ep, is estimated by the ratio of the 
measured peak stress and peak strain data. 

Traditionally, the root deformation has two stages, linear and 
non-linear or elastic and plastic deformation. The point that 
separates the two deformations is called the yield limit. It was 
commonly accepted that before the yield limit, the root 
deformation is linear or elastic and after the yield limit, the 
deformation is non-linear or plastic. However, data of our tensile 
tests show there is no clear point that separates linear or non-linear 
deformation. Instead, the dominated linear deformation is 
observed in early stage of tests, when the stress and strain are 
small. The weight of non-Hnear deformation gradually increases. 
At certain point, the non-linear deformation become significant 
and cannot be neglected. In this study, the term of the critical 
point is used to describe that point. Before the critical point, the 
non-linear deformation is very small when compared to the linear 
deformation (say less 5% of linear deformation). As the curve 
passes the critical point, the non-linear deformation becomes more 
and more significant. The critical point, Xp, is determined using 
Equation (2) and corresponding to the 40% of the tensile stress 
limit. 

Once above parameters are determined using experiment data, 
the proposed model is used to calculate the stress-strain relation. 
The error between the model calculated values and experiment 
data is estimated using the residual mean and variance and 
summarized in Table 2. The most average residual values are less 
than 0.1 and the most residual variance values are less than 0.01 
with only a few exceptions. Figure 3 presents eight figures of model 
calculated and experiment stress strain data. Seven of eight figures 
have a good match between model calculated values and 
experiment data, except the figure for Latrix gmelinii (diameter 
3.98 mm). Table 2 also shows that an increase in root diameter 
corresponds to a decrease in elastic modulus {Eg), while the same 
trend between the root diameter and the peak secant modulus {Ep) 
is not obvious. As results, the changing trend of E versus the root 
diameter is not clear. Experiment data with different diameters 
and species may be needed to further verify these observations. 

In summary, a constitutive model for the stress-strain relation- 
ship of single roots is proposed in this study. The model starts with 
a general form of linear and non-linear deformation and is defined 
with four boundary conditions based on typical stress-strain curves 
observed from the tensile tests of single roots. The definition of the 
model requires two physical characteristics of single root: the 
linear elastic modulus Eg and the peak secant modulus Ep. These 
two parameters can be obtained from the laboratory tensile 
strength tests. The major difference between the proposed model 
and most of previous stress-strain models [13; 16-1 7] is that the 
proposed model is based on mathematical expression of physical 
behavior of single root deformation and boundary conditions of a 
typical tensile test, instead of on experiment data only. The 
proposed model has a great potential to help the study of the 
mechanism of soil-reinforcement by root system. 
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Limitations 

The mechanism of soil-reinforcement by root systems is a 
complex issue and factors that may affect it include 1) 
environmental factors: such as species, direction and position of 
slope, soil properties, and root moisture content; 2) experimental 
factors: such as loading speed, loading time, length and diameter 
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Figure 3. Comparison between model calculated values and 

3.03 and 2.1 are root diameters. It can be seen that the theoretical ; 
doi:1 0.1 371 /journal.pone.0093066.g003 
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periment data of stress of selected testing. 1.1 1, 3.06, 2.15, 3.98, 1.06, 
!ss-strain curves basically fit the experimental curves well. 



of root, storage time of roots and root with or without bark; 3) 
micro-structural and chemical-composition factors: such as lignin 
and cellulose content, area ratio of wood fiber and phloem, and 
length-width ratio of wood fiber and tracheid. Our model 
developed in this study is based on laboratory root tensile tests. 
Although the proposed model has clear physical meanings and 
well defmed boundary conditions, it is limited by test data and 
conditions. Following is discussion about some limitations. 



3.1 Root Diameters and Tree Species 

As discussed above, our tensile tests are carried out in a range of 
root diameters (1.1 to 4.39 mm) of four tree species (four kinds of 
tall trees commonly found in Northern China) in one gauge 
(50 mm) used. The results presented in this paper are limited by 
those root diameters and tree species. The proposed model needs 
to be further verified with test data before it is apphed to different 
tree diameters and species and by different methods. 
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3.2 Loading Speed 

The loading speed used in our experiments was a constant speed 
at 1 0 mm/ min. Generally speaking, roots deform elastically before 
plastically. The root elastic deformation may not be fully 
developed if the loading speed is too high. Therefore, different 
loading speeds could affect the stress-strain relationship or the 
shape of curve. At extreme high loading speed, roots may suffer 
brittle fracture or necking in early stage of the test. 

3.3 Storage Time and Moisture Loss 

In our experiments, the maximum storage time was 30 days in 
order to avoid the moisture loss. Our data show that the moisture 
loss is insignificant within 30 days. The maximum moisture loss 

was ,3% in weight of roots for all roots and 1.4% for roots with 
diameters larger than 3 mm at the end of 30 days. The moisture 
loss could be significant after 7 weeks for Betula platjphylla, 1 1 weeks 
for Quercus mongolica, and 8 weeks for Ladx gmelinii and Pinus 
tabulaeformis. 

3.4 Witii or Without Barl< 

Another reason for the maximum storage time of 30 days is to 
prevent the deterioration of roots. It is assumed that the 
deterioration of roots will happen from the bark. Our test samples 
were with bark and no deterioration was observed in our test 
samples at the time of testing. The model developed in this paper 
is applicable for roots with bark. However, root samples with and 
without bark could make some differences in the stress-strain 
relationship, which subjects to further study. 

3.5 Boundary Conditions 

The four boundary conditions discussed in Section 2.3 describe 
a typical stress-strain relation of single roots, such as observed in 
our tests. The third boundary condition indicates the peak stress 
and strain (y= 1, when x — I) happen at the same time. The root 
breaks at the peak stress. This is true for almost all of test samples 
with a few exceptions. Some of test samples, although the number 
is small, do show the root deformation continued when the stress 
reached the peak value or necking. Figure 4 presents two stress- 
strain curves of such exceptions. It shows the peak stress does not 
correspond to the peak strain. The value of x=l does not 
represent the peak strain, instead it only represents the strain 
corresponding to the peak stress. The validity of the proposed 
model for this type of stress-strain relation is to be verified. 

Additionally, the elastic modulus. Eg, is determined by the 
secant modulus at 40 '/o of the tensile stress limit for all four tree 
species in this study. The selection of the secant modulus is based 
on previous study in which the elastic modulus ot Pinus tabulaeformis 
root system is 50% to 70% of the tensile strength limit. This study 
applies the same method to determine the elastic modulus of other 
three tree species. Whether the selection of 40% can be used for 
determination of the elastic modulus for different tree species used 
is out of the scope of this study and deserves further study. 

Even though with some limitations, the proposed model is 
promising, since it's deduced from mathematical procedures and 
verified by a large amount of test data. The application of the 
model is simple. It is one of first constitutive models of single root 
resistance to tensile stress with both physical and mathematical 
meaning. Systematic and detailed studies on above issues are 
subjects of further research. 

Conclusions 

Although the ultimate aim of our research is to study the role of 
tree roots in the slope stabilization, the specific aim of this paper is 
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Figure 4. Examples of roots that broke after the stress reached 
the peak value. 
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to build a constitutive model of single root system with a clear 
physical meaning through the rigorous mathematical derivation. 
The development of the constitutive model and the study of single 
root physical response to tensik- stress are fundamentally 
significant to the understanding of soil-reinforcement by root 
systems. Since vegetation is one of major factors that affect slope 
stability, a common question is how exactly vegetation (especially 
the mechanical property of roots) affects the soil stabilization. The 
proposed model of single root resistance to tensile stress and the 
study reported in this paper provide a theoretical basis for the 
study of the slope stabilization by vegetation. 

The mechanism of soil-reinforcement by root system is 
investigated by studying single roots of four common tall tree 
species from Northern China. The data of laboratory tensile tests 
is used to construct the root tensile constitutive model. The 
following is a summary of this study: 

5) A constitutive model of single root system resistance to tensile 
stress is proposed. 

6) The proposed model is mathematically defined by four boundary 
conditions observed in a typical tensile test of single root. 

7) The term of the critical point is introduced. Before the critical 
point, the linear deformation is dominated and after the critical 

point, the non-linear deformation becomes significant. 

8) The physical parameters of single root systems, linear elastic 
modulus [Egj and peak secant modulus {Epj, can be determined by 

laboratory tensile tests. 

9) Large amount of experiment data from the tensile tests of four 
common tall trees' root systems are used for the model verification. 
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The residual mean and variance between the model calculated 
values and experiment data indicates the proposed model can be 
used for the study of stress-strain relationships of single roots with 
good accuracy. 

10) The proposed constitutive model of single roots provides 
theoretical basis for the research of the mechanisms of soil- 
reinforcement by root system and makes the quantitative analysis 
of the root tensile strength feasible. 
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